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Abstract 

Using data taken with the CLEO-c detector at the Cornell Electron Storage Ring, we have 
investigated the direct photon momentum spectrum in the decay J/ip(lS) — > 'ygg, via the "tagged" 
process: e + e~ — > ip(2S); ^(2S)— > J/tp ir + ir~; J/ip — > 7+X. Including contributions from two-body 
radiative decay processes, we find the ratio of the inclusive direct photon branching fraction to 
that of the dominant three-gluon branching fraction (Ry = B(ggj)/B(ggg)) to be i? 7 = 0.137 ± 
0.001 ± 0.016 ± 0.004, where the errors shown are statistical, systematic, and the model-dependent 
uncertainty related to the extrapolation to zero photon energy. The shape of the scaled photon 
energy spectrum in J/tp — > gg^ is observed to be very similar to that of T — > ggj. The Ry value 
obtained is roughly consistent with that expected by a simple quark-charge scaling (Ry ~ (ic/lb) 2 ) 
of the value determined at the T(1S), but somewhat higher than the value expected from the 
running of the strong coupling constant. 

PACS numbers: 13.20.Gd,13.20.-v,13.40.Hq 
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I. INTRODUCTION 



According to the Okubo-Zweig-Iizuka (OZI) rule, the preferred decay mode for char- 
monium would be through the production of a DD meson pair. For resonances below 
the V(3770) however, this is not energetically possible. Thus the decay of the J/ip(lS) 
meson must proceed through OZI-suppressed channels. The three lowest-order decay 
modes of the J/ip meson are the three-gluon (ggg), virtual photon (vacuum polariza- 
tion) decays J/ip — > l + l~ and J/ip — > qq with a branching fraction given in terms of 
IZqcd = (e + e~ — > qq)/(e + e~ — > as ~ TZqcd x B(J/iIj IT), and two gluon plus 

single photon (ggj) modes. For the ip(2S) and -?/>(3770) resonances, direct radiative transi- 
tions, both electromagnetic and hadronic, as well as decays to open charm (for the t/>(3770)), 
compete with these annihilation modes and therefore reduce the gg^ branching fraction. 



A. Inclusive Total Rate 



Since Y ggg oc and T ggi oc a 2 a em , the ratio of the branching fractions for the ggg and 
<7<?7 decay modes for heavy quarkonia (also equal to the ratio of experimentally measured 
events N ggi and N ggg , respectively) is expected to follow [l[: 

Ry = = ^ = T Vc — i 1 + (2-2 ± 0.6)a s M (1) 
a ggg ^ggg a $ 

In this expression, the charm quark charge q c = 2/3. Alternately, one can normalize to 
the well-measured dimuon channel and cancel the electromagnetic vertex: B ggi /B^ oc a 2 s . 
In either case, one must define the momentum scale {Q 2 ) appropriate for this process. 
Although the value Q 2 ~ M 2 seems natural, the original prescription of Brodsky, Lepage 
and McKenzie ("BLM" Q) gave Q = 0.157M T( i S) for (the less-relativistic) T(lS)-> ggj. 
Alternative prescriptions for the appropriate value of Q 2 have also been suggested (2|. 



B. Energy and angular spectrum shapes 

Calculations of the direct photon energy spectrum were originally based on decays of 
orthopositronium into three photons, leading to the expectation that the J/ip direct pho- 
ton energy spectrum should rise linearly with z 1 (= 2E 1 /Mj/^) to the kinematic maximum 
(Zj — > 1); phase space considerations lead to a slight enhancement at z 7 =l jl]. The angular 
distribution for the decay of a polarized vector into three massless vectors is, in principle, 
directly calculable. Thus, for direct radiative decays J/if> — > gg^f, theory prescribes the 
correlation of z 7 with photon polar angle cos# 7 , defined relative to the beam axis. Koller 
and Walsh considered the angular spectrum in detail Q, demonstrating that, if the parent 
is polarized along the beam axis, then, as the energy of the most energetic primary (photon 
or gluon) in J ftp — * jgg or J ftp —>■ ggg approaches the beam energy, the event axis tends 
to increasingly align with the beam axis: z 1 — > 1 corresponds to a(z 7 ) —>■ 1 for an angular 
distribution specified as dN/d cos # 7 ~ 1 + a(z 7 )cos 2 # 7 . We note that, according to the 
Koller-Walsh prescription, the value of a(z 7 ) for intermediate values, where most of the 
events occur, is relatively small (0.2). Only for z > 0.9 is the forward peaking of the photon 
angular distribution noticeable. 
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Previous analyses of the direct photon spectrum in heavy quarkonium decay selected a 
fiducial angular region and integrated over cos# 7 . In this analysis, we will take advantage 
of the expected correlations between cos# 7 and to improve the statistical precision of 
the extracted branching fraction. There is nevertheless still some model dependence in the 
extrapolation down to z 7 — > 0. 



C. Previous Work 

Garcia and Soto (GS) [sj have performed the most recent calculation of the expected 
direct photon spectrum in J / ip decays, using an approach similar to that applied by the same 
authors for the case of T(1S)— > gg^ |6]. They model the endpoint region by combining Non- 
Relativistic QCD (NRQCD) with Soft Collinear- Effective Theory (SCET), which facilitates 
calculation of the spectrum of the collinear gluons which occur as z 7 — *■ 1. Both color-octet 
and color-singlet contributions must be explicitly calculated and summed. The calculations 
are very sensitive to the handling of the octet contribution, and limit the momentum interval 
over which the theory is considered 'reliable' to 0.4 < z 7 < 0.7. At low energies (defined as 
z 7 < 0.45 for the case of the T), the so-called "fragmentation" photon component, due to 
photon radiation from final-state quarks, dominates. 

Although inclusive radiative decays have received considerable experimental attention in 
the case of bb 0, @, 0, H3, U, 13], the cc system has had only one prior measurement, by the 



MARK-II Collaboration in 1981 [13(. The MARK-II analysis, which utilized a calorimeter 



with resolution /Je/E ~ 0.12/y E(GeV) resulted in a measurement for the inclusive partial 
branching fraction B(J/ip — *> r yX)={4.1 ± 0.8)% limited to the range z 7 >0.6. Although the 
authors do not explicitly quote a value for i? 7 in their original reference, we can estimate an 
implied i? 7 value assuming that z 7 > 0.6 constitues 45% of the total spectrum (over the full 
cos# 7 range, based on results obtained for bottomonium), and using T^qcd =31 GeV =2.1 and 
8^=0.059, such that B(ggg) + B(ggj)=0.69, yielding R, = B{ ggi )/B{ggg) ~ (14.6±2.8)%. 
The MARK-II direct photon spectrum peaked at z 7 ~ 0.6, inconsistent with expectations 
based on orthopositronium decay, but consistent with later bottomonium spectra. 



II. DETECTOR AND EVENT SELECTION 

The CLEO-c detector is essentially identical to the previous CLEO III detector, with 
the exception of a modified innermost charged particle tracking system. Elements of the 



detector, as well as performance characteristics, are described in detail elsewhere [14J, [15|, [16 
Over the kinematic regime of interest to this analysis, the electromagnetic shower energy 
resolution is approximately 2%. The tracking system, RICH particle identification system, 
and electromagnetic calorimeter are all contained within a 1 Tesla superconducting solenoid. 

In the absence of dedicated J ftp data collected at CLEO-c, we use the cascade decay 
chain ip(2S) — > J/ip7i + n~; J/ip — > ggj. Our data sample corresponds to approximately 
27xl0 6 ip{2S) decays[I3] collected with the CLEO-c detector (our "primary" data sample, 
divided into two sub-samples of data taken approximately three years apart); a much smaller 
("secondary") sample of slightly more than 10 6 ip(2S) decays collected with the CLEO-III 
detector is used for cross-checks. To ensure maximal efficiency, minimal event selection 
requirements are imposed on our candidate direct photon sample - we require only that 
candidate events have at least two high-quality charged tracks (the two transition charged 
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pion candidates) and no identified lepton charged tracks. Our lepton veto is effective in 
suppressing contamination from J/if> — > . Monte Carlo simulations indicate that the 
trigger efficiency for events having two transition charged pions and a fiducially-contained 
direct photon with z 7 >0.3 is >99%. 

III. ANALYSIS 

To obtain Rj, we must determine separately the number of direct photon events and the 
number of three gluon events. By using the decay chain if>(2S) — > J/ip7r + 7r~; J /if} — > gg'y, we 
circumvent initial state radiation backgrounds in our photon sample. The dominant back- 
ground for z 1 <0.6 is primarily from tt° — > 77 and 77 — *> 77. Photon selection requirements 
are essentially the same as applied for our study of T — > gg^f [?J . Namely, we require showers 
detected in the electromagnetic calorimeter with energy deposition characteristics consistent 
with those expected for true photons, and which are well-isolated from both charged tracks 
as well as other showers. Since the 7r° signal to noise is high, we will, in contrast to our 
previous analysis, require that a candidate high-energy photon not combine with another 
high-quality photon to give an invariant mass within ±6 MeV from the nominal tt° mass 
m°. However, we do not impose an rj veto given the worse signal to noise and the greater 
likelihood of incorrectly vetoing a true direct photon. The scatter plot of the raw candidate 
photon energy vs. "shifted" dipion recoil mass (ECM — M recoi i) is shown in Figured) 
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FIG. 1: Center-of-mass energy minus dipion recoil mass (horizontal) vs. photon energy (vertical), 
V>(2S) data. 

This spectrum includes contributions from each of the three main decay modes of the J/if>: 
three-gluon, vacuum polarization, and direct photon. Photons from nonresonant processes 
below the resonance, e + e~ — > qq^), also contribute to this spectrum. However, since the 
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charmonium peak cross-section is a factor of 50 times larger than the continuum cross- 
section, and we perform a dipion-sideband subtraction to obtain the "tagged" direct photon 
spectrum from the J/ip, we will, in what follows, ignore this continuum contribution. For 
the sideband subtraction, the signal region is defined as a dipion recoil region within ±10 
MeV/c 2 of the canonical J/ip mass; sidebands are defined as the regions 10-40 MeV/c 2 from 
the canonical J/ip mass. 

Knowing the dipion four- vector allows us to work in the rest frame of the J / ip itself; in 
what follows, unless otherwise indicated, the energy spectra presented correspond to this 
case. 



IV. BACKGROUND PHOTON SUBTRACTION AND SIGNAL ESTIMATE 

In the previous bottomonium analysis, two parameterizations of the background were used 
- one was based on the "pseudo-photon" technique which has been used in three previous 
CLEO analyses @, I, H^ |. as well as the original MARK-II analysis, and the other used a 
simple exponential parametrization of the background under the direct photon signal. The 
latter suffers from integrating over the correlations between cos# 7 and z 1 and is therefore 
not used in our current analysis. To model the production of ir° and rj daughter photons, 
redundant estimators were employed in this analysis, which we detail below. Unlike the 
bottomonium analysis, we do not simulate uj — > 7r°7, and v( — > j(p,u,j) contributions. 
JETSET M indicates that these should be smaller than for the T. Numerically, the fraction 



of all 2 7 >0.4 photons having uj{rf) parentage is 2.2%(0.8%) according to the LUND event 
generator, somewhat below our typical systematic errors. The non-photonic contribution to 
our final candidate shower sample, due almost exclusively to Kl and n interactions in the 
calorimeter, is estimated from Monte Carlo simulations to be 1.5%. We note that, since the 
background shape is fixed during signal extraction, but the normalization allowed to float, 
such non-direct photon contamination is largely absorbed into the eventual background 
estimate. 



A. "Pseudo-photons" 

As a first estimate of the non-direct photon background, we took advantage of the ex- 
pected similar kinematic distributions between charged and neutral pions, as dictated by 
isospin symmetry. Although isospin will break down both when there are decay processes 
which are not isospin-symmetric in their final states {J/ip — > 777, r\ — > 7r°7r 7r°) or when the 
available fragmentation phase space is comparable to M n , in the intermediate-energy regime 
we expect isospin to be reliable. Over the kinematic regime relevant for this analysis, Monte 
Carlo studies indicate consistency (to within ~5%) with the naive expectation that there 
should be half as many neutral pions as charged pions, with similar momentum-dependent 
angular distributions. However, unlike the case for the T, the region z 1 — >1 has large contri- 
butions from two body radiative decays. Specifically, the 7r° : ratio grows in this regime 
due to decays such as J/ip — > 777 (rj — > 37r°) and J/ip — > 777'. 

A pseudo-photon spectrum (d 2 N/ (dz 7 d cos 9 7 )) is constructed using charged tracks with 
particle identification information consistent with pions to model the spectra expected for 
7r°'s and 77's. In both cases, we use the Monte Carlo prescribed tt° : ^ or rj : 7r ± ratios, 
taking into account the variation in these ratios with momentum. Momentum-dependent 
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corrections are also applied to account for non-pion charged kaon, proton, and lepton fakes 
in our sample, as well as the finite charged track-finding efficiency. Each of these last two 
corrections are of order 5% and tend to offset each other. We invoke isospin in making the 
assumption that the momentum-dependent angular distributions of neutral pions follows 
that of charged pions. Our Monte Carlo generator indicates that the momentum-dependent 
angular distribution for 77's is also similar to charged tracks. We simulate the two-body 
decays tt° — > 77 and 77 — > 77 in the rest frame of the candidate tt° or 77 parent and boost 
the daughter pseudo-photons into the lab frame according to the (tt° or 77) momentum. The 
direct-photon finding efficiency e 7 (.E 7 ) is applied to each daughter pseudo-photon to deter- 
mine the likelihood that that photon will populate our candidate direct photon spectrum. 
From Monte Carlo simulations, e 7 ~ 0.85 over the kinematic and geometric fiducial interval 
defined in this measurement. Finally, "found" pseudo-photons are smeared in energy and 
angle by the known resolutions. 

To check our procedure, we have compared the data 77 invariant mass plot with the 
pseudo-photon 77 invariant mass plot. To enhance statistics, we use all photons with 
E 1 >0.2 GeV {z 1 > 0.13). At such relatively low photon energies, the number of accepted 
showers not having tt° or 77 parentage is considerable, so we must also add to our M 77 
spectrum combinations of "found" n° or rj daughter pseudo-photons with "excess photons" , 
using a Monte Carlo-based momentum- dependent factor. We obtain a level of agreement 
(better than 10%) consistent with our previous T(1S) analysis |7j. The comparison between 
our pseudo-photons and data is shown in Fig. [2j 
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FIG. 2: Two-photon invariant mass combinations in data (crosses) vs. pseudo-photon simulation. 
Shaded histogram indicates pseudo-photon simulation after adding combinations of simulated ir° 
and 77 daughters with "excess" showers in data. Normalization is absolute. 
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B. Background Estimate from Monte Carlo Simulations 

Second, we use the Monte Carlo simulation of generic J/ ip decays, based on the JETSET 
7.4 event generator, to provide an estimate of the background to the non-direct photon 
signal, including all sources. This estimate implicitly includes all the corrections (photon 
efficiency, tracking efficiency and fake rates, tt° : 7r =l= ratio, hadronic showers, etc.) which must 
be explicitly evaluated in the previous approach. No additional corrections are therefore 
applied in this case. 



C. Background using XcO decays 



A third estimate of the background is obtained from decays of the Xco(3415), which is 
produced via the transition: ip(2S)-^ Xcol with the emission of a photon of energy 270 
MeV. Since the Xco is relatively wide (r ~10 MeV), it has a relatively small radiative 
decay branching fraction to the J /if) (23=1.32 ± 0.11%) and therefore dominantly decays 
via two-gluon intermediate states. To the extent that two-gluon fragmentation is similar to 
three-gluon fragmentation 191 ]. we can therefore use the data photon background produced 
in association with an observed if) (2S)—> Xcol transition photon candidate to estimate the 
non-direct photon background to the if) — *> gg'y photon energy spectrum. Comparison of 
the charged track spectra for sideband-subtracted Xco vs. J /if) decays indicates that the 
kinematics of the former two-gluon decays are similar to the latter three-gluon decays in 
this case (Figure [3]). 




-| Q2 I I ! ! I ! ! ! ! I ! ! ! ! I 

0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 

Charged track momentum (GeV/c) 

FIG. 3: Comparison of XcO —> ^ (shaded) vs. J /if) — > 7r (dashed) inclusive charged track spectra. 

To further suppress any possible Xco J /''Pi cascade contamination, we veto events 
which contain a high-quality photon candidate of energy 322±20 MeV. A signal region is 
defined around the 240 — 290 MeV transition photon energy range and the photon energy 
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spectrum in coincidence with the candidate transition photon (after sideband subtraction, 
with sidebands taken an additional 25 MeV on each side of the signal region) is then used 
as our last background-estimator for our final fits. 

D. Polarization of the parent J/ip 

In principle, the dipion transition can be either S- or D-wave, as allowed by parity conser- 
vation. The BES Collaboration have studied the angular distributions for this process [2(| 
and find a best-fit value for the D-wave to S-wave amplitude of 0.18±0.04. If the decay is 
all D-wave, then the 1 + a cos 2 # 7 distribution expected for two-body 7+pseudoscalar decays 
softens to a ~ 0.07 compared to a = 1 for S-wave. (We use the symbol a to designate 
the angular distribution of the radiative daughter photon specifically in two-body radiative 
decays.) Similarly, there is some uncertainty in the angular distribution of the direct photon 
signal itself, characterized by the inclusive spectral parameter a(^ 7 ). To accommodate this, 
we have done fits varying values of both the two-body angular coefficient ct(z 7 ) and the 
inclusive direct photon angular coefficient a, and include the difference among them as a 
systematic error. 

E. Fits and Signal Extraction 

After imposing our event selection and photon selection criteria, we are left with the 
two-dimensional candidate direct photon scaled energy vs. polar angle distribution. We 
perform two-dimensional fits comprised of the following components: a) the background, 
which is modeled either using the pseudo-photon, Monte Carlo-based, or x c0 -based back- 
grounds described above, b) three two-body components of the direct photon signal: 777 
and 777' and a wide resonance which corresponds to 77/(1440), with shapes determined from 
Monte Carlo simulation, and c) a smooth signal component which has a shape in photon 
energy taken from our previous Upsilon decay measurement and an angular distribu- 
tion based on the Koller- Walsh prescription [4J. Ideally, we could avoid having to include 
a signal component. In such a case, the background subtraction would directly determine 
the true underlying signal. However, this can only be done if the background can be abso- 
lutely normalized with very high precision (much better than our ~10% overall background 
normalization error). Unfortunately, the statistics of the fit are largely driven by the low-z 7 
region, where the systematic uncertainty on the 7r° background is largest. Without inclusion 
of a signal component, the background normalization would increase to saturate the low-27 
energy regime. 

A comparison of the one-dimensional projections of the background dN/dz^ spectra is 
shown in Figure HI 

F. Validation of fit procedure using Monte Carlo simulations 

For the Monte Carlo simulations, we can fit the signal plot with a combination of the 
tagged Monte Carlo simulation signal plus either the pseudo-photon or the Xco-based back- 
ground, then check the signal normalization against the known number of signal direct 
photons in the Monte Carlo simulations. The Xco background estimator results in a signal 
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FIG. 4: Comparison of background photon {{l/N)dN/ dz~) spectra for three subtraction schemes 
considered. Plus signs: pseudo-photon background; histogram: XcO background; circles: Monte 
Carlo simulations of background. 

yield only 3% larger than the known number of Monte Carlo direct photons in the plot, while 
the pseudo-photon background estimator underestimates the signal yield by ~15%. We note 
that the agreement between the signal yields obtained from these two backgrounds in data is 
typically within 3%. As discussed later, we nevertheless add an additional systematic error 
(6%) to reflect this discrepancy observed in simulation. 

V. RESULTS 

Figs. O [6] and [7] show sample fits over the kinematic region z 7 >0.3 and | cos# 7 | <0.8, 
based on the pseudo-photon model of the background, the Monte Carlo-based model of the 
background, and the Xco-based model of the background, respectively. 

Positive residuals (data in one bin exceeds sum of fit contributions) are shown in black; 
negative residuals (data in one bin is smaller than sum of fit contributions) are shown 
in white. Individual projections, onto the scaled photon energy and photon polar angle 
axes, for the three background models separately, are presented in Figure [SJ An overlay 
of the background-subtracted spectra, for the three background models employed, is shown 
in Figure [9] for our primary dataset. We observe reasonable agreement between the three 
spectra over most of the kinematic regime considered. 

Excluding electromagnetic transitions to other charmonium states, two-body radiative 
exclusive modes should be included in our total gg^ yield. Decays into narrow r\ mesons 
have large enough branching fractions (#(777) = (9.8 ± 1.0) x 10~ 4 and B{p(r]') = (4.71 ± 
0.27) x 10~ 3 ) so that they are clearly visible in the one-dimensional projection of the signal 
spectrum. Rather than fixing these contributions, we have allowed them to float and use the 
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FIG. 5: Fit over z 7 > 0.3, | cos# 7 | <0.8 kinematic region using pseudo-photon background. Shown 
are normalized residuals, in units of statistical error per bin. Black: data>fit; White: data<fit. 
Excess at large values of cos # 7 attributed to QED processes producing charged leptonic tracks at 
large dip angles, which are (incorrectly) used as input to pseudo-photon generator. This results in 
excess pseudo-photons at high values of cos6> 7 . 




FIG. 6: Fit over z 7 > 0.3, | cos# 7 | <0. 8 kinematic region, MC background. Black: data>fit; White: 
data<fit. Note apparent presence of extra two-body component in data, at z 7 « 0.5 — 0.55. 
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FIG. 7: Fit over z 7 > 0.3, | cos 9j\ <0.8 kinematic region, XcO background. Black: data>fit; White: 
data<fit. We again note apparent excess of data over background in high cos # 7 region around 



0.4 < z 7 < 0.5. 



fitted area, corrected by efficiency, as a check of the overall procedure. We obtain estimates 
for the two-body branching fractions of (8.1 ± 0.6) x 1CT 4 and (4.98 ± 0.08) x 10~ 3 for our 
primary data sample and (8.8 ± 2.1) x 10~ 4 and (5.4 ± 0.4) x 10 -3 for our lower statistics, 
secondary data sample, where the errors presented are statistical only. Systematic errors for 
this coarse cross-check are likely to be at least as large. 



VI. EXTRACTION OF # 7 

To determine _R 7 , we calculate the ratio of true N ggi and N ggg events (Eq. (1)). The 
former quantity is the number of observed direct photon candidates, corrected by the photon- 
finding efficiency and the kinematic acceptance. The latter quantity is the total number of 
^(2S)— > J/ipTT + TT^ events estimated from our inclusive dipion recoil mass spectrum, minus 
contributions from J/ip — > gg^ (implicitly including all two-body radiative components), 
J /ip —> / + /~, and J/ip — > qq. The two-track trigger efficiency is >99% for all these processes. 
The hadronic event reconstruction efficiency is also >99%; due to our lepton veto, the 
acceptance for dileptonic decays is less than 5%. 

Table [T] shows the range of values obtained for two sub-samples of the primary data, taken 
approximately three years apart (presented as "sub-sample-l/sub-sample-2") with varying 
definitions of the signal region. 

VII. SYSTEMATIC ERRORS 

We identify and estimate systematic errors in our i? 7 determination as follows: 
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FIG. 8: Background-subtracted signal yield projections onto z 7 and cos# 7 axes, for the three 
photon background estimators used in this analysis. 

1. The uncertainty in the S:D admixture of the dipion transition, in principle, affects the 
angular distribution (~ 1 + a cos 2 9 y ) and therefore the acceptance for the two-body 
radiative component in our fits, as well as the angular distribution, and acceptance for 
the primary direct photon signal component (~ 1 + a(z 7 ) cos 2 # 7 ). Although both S- 
wave and the dominant allowed D-wave transition amplitudes leave the daughter J / ip 
polarized along the beam axis, we nevertheless allow for possible contributions due to 
D-wave amplitudes resulting in the daughter J/ip polarized transverse to the beam 
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FIG. 9: z 7 > 0.35, | cos 9y\ <0.8, background-subtracted direct photon energy spectra, using three 
different background subtraction schemes. Also overlaid (shaded histogram) is the experimental 
spectrum for the T spectrum. 



TABLE I: Summary of values of i? 7 using varying kinematic regions for fits, different datasets and 
background subtraction schemes. First number is primary sub-sample 1, second is primary sub- 
sample 2. Numbers presented are simple averages of the values obtained using the three different 
background estimators. The number in parenthesis following each average corresponds to the rms of 
the difference among the three background estimators and is therefore indicative of the (dominant) 
signal estimation systematic error. By comparison, statistical errors are typically of order 1% for 
sub-sample 1 and 0.2% for sub-sample 2. Our final quoted i? 7 result is a weighted average of the 
presented values. 

| cos6> 7 | <0.7 |cos# 7 |<0.8 
z 7 > 0.30 0.1367(88)/0.1398(78) 0.1338(97)/0.1371(85) 
z 7 > 0.45 0.1320(110)/0.1362(94) 0.1332(89)/0.1403(98) 



axis. Varying a between 0.7 and 1.0 for the two-body modes results in a 1% change 
in the extracted direct photon yield; varying a(z 7 ) between the value prescribed by 
Roller- Walsh and a=0 for all photon momenta results in a 3% lower value for the 
extracted direct photon yield. We attribute this to the larger saturation of the signal 
region by background which results when the signal photon angular distribution is 
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taken to be flat in angle. We assume a systematic uncertainty of 3% due to this 
source. 

2. The uncertainty in the contribution to the signal due to non-photon showers, based 
on Monte Carlo modeling of signal and background decays, is estimated to be 1.5%. 

3. Uncertainty in the number of 3-gluon events is obtained by subtracting from the 
total number of dipion tags the number of J/ip — > 7* — > qq events, the number of 
signal J/ip — > 7* — > gg^f events, and the number of J/ip — > 7* — > dileptonic 
decay events which pass our cuts. The statistical error on the branching fraction for 
J/ip — > 7* — » qq is very small (£> = 13.50 ±0.30%) |2l| . as is the error on the dileptonic 
branching fraction (5.94 ± 0.06%). The fraction of dileptonic events which pass our 
cuts is also small (<5%), as is the statistical error on the number of gg'j events in our 
sample. The total systematic error on our calculated ggg yield due to the noia-ggg 
subtractions is largely due to the uncertainty in i? 7 and determined to be ~2%. 

4. The trigger efficiency systematic error in the ratio is <1%. 

5. Background normalization and background shape uncertainty are evaluated by exam- 
ining the agreement between the direct photon yield obtained using the three different 
background estimators. We point out that these three techniques sample very dif- 
ferent methods of background-estimation. The Xco subtraction background estimate, 
e.g., is insensitive to the uncertainty in the overall photon-finding efficiency Given the 
observed agreement across momentum, we infer that the pseudo-photon technique is 
least likely to be sensitive to ir°/r] modeling uncertainties. 

However, we observe that the fits follow a generally consistent pattern. Although the 
data-driven fits (pseudo-photon and using the Xco background) are generally consistent 
with each other at the 3% level, the average of the data-driven fits is consistently lower 
(by ~15%) than the Monte Carlo-background subtracted spectra. The overall rms of 
the signals obtained using the three background-estimators is 7% (see Table [T|) and we 
assign this value to the corresponding systematic error. 

6. The uncertaingy in the absolute photon-finding efficiency is estimated at 2%. 

7. Sensitivity to the selection of signal and sideband regions in the dipion mass spectrum 
is estimated by increasing the nominal 'signal' recoil mass interval by 25% and de- 
creasing the nominal 'sideband' recoil mass interval by 25%, indicating a systematic 
error <2%. 

8. The difference in our calculated value of i? 7 between imposing vs. not imposing the 
7T° veto is found to be about 3%. 

9. For the pseudo-photon subtraction only, the sensitivity to the assumed 7r° : tt ratio 
was estimated by comparing the results based on the Monte Carlo-prescribed ratio vs. 
a constant value of 0.5. This results in a variation of 3% in R 1 . 

10. Possible continuum QED contamination should be subtracted out via the dipion side- 
band subtraction, although we do rely on Monte Carlo simulations to quantify the 
background from processes such as ^>(2S)— >• J/tpn+n^; J/ip — > Our results 

with very strict QED suppression vs. no QED suppression vary by 1%. We conserva- 
tively assign a 1% systematic error due to our uncertainty in this background. 
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11. As described previously, we have compared the signal yield with the "true" signal yield 
using a Monte Carlo-only study, in which the number of simulated signal photons are 
known. Unfortunately, the JETSET 7.4 Monte Carlo simulated spectrum is entirely 
two-body and quasi-two-body, and does not reproduce data well (Figure [TU]). As 
outlined previously, we find that our average extracted signal yield is smaller than the 
true signal magnitude in Monte Carlo simulations by 6%, and conservatively (since this 
error likely is somewhat redundant with the systematic error assessed by the spread 
in i? 7 values obtained using the three different subtraction schemes) include this as an 
additional systematic error. 
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FIG. 10: Monte Carlo simulations of direct photon J/tp radiative decay photon spectrum. 



A separate, additional systematic error must be included to account for the uncertainty 
in the extrapolation to 2 7 =0. For this, we compare the values obtained assuming a linear 
extrapolation from z 7 =0.3 to z 7 =0 vs. an extrapolation based on the shape of the spectrum 
observed in the case of the T. The difference between the yields for these two extrapolations 
is <3%. 

Table Hn compares the results of this analysis with those obtained by previous experiments. 
Although the statistics are poorer, the older CLEO-III data (our cross-check sample) gives 
results which are consistent with the CLEO-c results (0.132 ± 0.008 ± 0.013, where the first 
error is statistical and the second represents the spread in the measured values obtained 
using the three different background subtraction schemes). 



VIII. IMPLICATIONS FOR a s 

Although the large relativistic corrections may render such an estimate unreliable, we can, 
nevertheless, calculate the value of a* implied by our i? 7 measurement, as shown in Figure 
[TTi Voloshin, in his recent review [22j, estimates an expected branching fraction B(J/ip — > 
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Experiment 



MARK-II [12] 0.041 ± 0.008 (z 7 > 0.6 only) 
0.146 ± 0.028 (all z 7 , estimated) 
This measurement R~ = 0.137 ± 0.001 ± 0.016 ± 0.004 



TABLE II: Comparison with previous experiment. MARK-II errors are total. CLEO-c errors are 
statistical, systematic and the uncertainty in the extrapolation to zero direct photon momentum. 
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FIG. 11: Horizontal band indicates R~ central value and experimental uncertainties; curve cor- 
responds to predicted a s band (including theoretical errors), using the BLM prescription for 
Q 2 = 0.157M 9? . 



gg r y)=6.7%, using a s (m c )=0.19 and the known value of T ee (J/ip). We can translate our 
value of Rj into B by correcting for the non-ggg decays of the J /if) (65.5%), to obtain 
B(J/if> — > ggj) = 9.2 ± 1.0%, considerably higher than Voloshin's estimate. We note that 
the earlier MARKII result, extrapolated to the full kinematic regime is also somewhat larger 
than Voloshin's estimate. 



IX. SUMMARY 

We have extracted the direct photon energy spectrum from J /if) decays based on a two- 
dimensional fit procedure. Normalized to the dominant 3-gluon mode of the J /if), and 
including two-body radiative decays, we obtain i? 7 = 0.137± 0.001 ±0.016 ±0.004. Although 
consistent with the one previous measurement, our direct photon yield is somewhat higher 
than that expected by a simple extrapolation from results at the T(1S) (-R 7 = 2.78 ±0.08%, 
averaged over all previous measurements), taking into account the variation in a s {Q 2 ). 
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